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ABSTRACT:  Development  of  repeatable  and  quantita¬ 
tive  tools  are  necessary  for  determining  the  abundance 
and  distribution  of  different  types  of  benthic  habitats,  de¬ 
tecting  changes  to  these  ecosystems,  and  determining 
their  role  in  the  global  carbon  cycle.  Here  we  used  ocean 
color  remote  sensing  techniques  to  map  different  major 
groups  of  primary  producers  and  estimate  net  primary 
productivity  (NPP)  across  Great  Bahama  Bank  (GBB). 
Field  investigations  on  the  northern  portion  of  the  GBB  in 
2004  revealed  3  dominant  types  of  benthic  primary  pro¬ 
ducers:  seagrass,  benthic  macroalgae,  and  microalgae  at¬ 
tached  to  sediment.  Laboratory  measurements  of  NPP 
ranged  from  barely  net  autotrophic  for  grapestone  sedi¬ 
ment  with  thin  microalgal  biofilm  to  highly  productive  for 
dense  accumulations  of  brown  macroalgae.  A  logarith¬ 
mic  relationship  between  NPP  and  green  seafloor  re¬ 
flectance  described  the  general  trend  in  NPP  across 
various  benthic  constituents.  Using  a  radiative  transfer- 
based  approach,  satellite-derived  estimates  of  NPP  for 
the  region  totaled  -2  x  1013  gC  yr_1  across  the  GBB.  The 
prevailing  benthic  habitat  was  mapped  as  sediment  with 
little  to  no  microalgal  biofilm.  Moderate  to  dense  sea¬ 
grass  meadows  of  Thalassia  testudinum  were  the  domi¬ 
nant  primary  producers  and  contributed  over  80%  of 
NPP  in  the  region.  If  the  vast  majority  of  seagrass  leaves 
decompose  in  the  primarily  carbonate  sediments,  car¬ 
bonate  dissolution  processes  associated  with  this  decom¬ 
position  may  result  in  sequestration  of  seagrass  above- 
and  below-ground  carbon  into  the  bicarbonate  pool  (2.4 
x  1013  gC  yr_1),  where  it  has  a  residence  time  on  the  order 
of  tens  of  thousands  of  years. 
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Pseudo-true  color  image  from  the  MODIS  Aqua  sensor  at 
250  m  resolution  shows  the  multihued  reflections  from 
shallow  waters  in  the  southern  Bahamas. 

Image:  C.  Buonassissi  and 
NASA  Ocean  Biology  Processing  Group 


INTRODUCTION 

In  1918,  Johannes  Petersen  published  the  first  spa¬ 
tial  estimates  of  seagrass  production  along  a  7000  km2 
basin  off  the  coast  of  Denmark  and  reported  that  it 
was  more  than  4  times  the  quantity  of  hay  annually 
produced  by  Danish  fields  and  meadows  (Petersen 
1918).  Despite  this  early  effort,  global  assessments  of 
productivity  from  seagrass  and  other  benthic  produc¬ 
ers  remain  elusive  and  their  contributions  to  regional 
and  global  biogeochemical  cycles  remain  poorly 
quantified  (Smith  1981,  Short  &  Wyllie-Echeverria 
1996). 
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Benthic  producers  can  have  relevance  to  climate 
studies  both  in  terms  of  the  amount  of  productivity 
and  the  fate  of  that  productivity.  For  example,  sea- 
grasses  may  cover  roughly  6  x  105  km2  of  the  global 
seafloor  and  their  net  primary  production  (NPP)  is 
estimated  to  be  0.5  PgC  yr_1  (Duarte  &  Cebrian  1996, 
Mateo  et  al.  2006).  This  reflects  roughly  1  %  of  NPP  of 
the  ocean’s  phytoplankton  (Longhurst  et  al.  1995). 
The  role  benthic  producers  play  in  both  the  local  ecol¬ 
ogy  and  biogeochemical  cycles,  however,  also  depends 
on  how  much  production  is  left  unremineralized  by 
herbivores  and  decomposers  to  accumulate  in  the 
seabed  (Short  1980,  Mateo  et  al.  2006).  Along  the 
Exumas  portion  of  the  Great  Bahama  Bank  (GBB), 
highly  productive  mats  of  negatively  buoyant  benthic 
macroalgae  may  be  transported  into  the  deep  Tongue 
of  the  Ocean  in  pulsed  exports  exceeding  7  x  1010  g  of 
carbon  (Dierssen  et  al.  2009b).  Such  events  are  equiv¬ 
alent  to  the  daily  carbon  flux  of  phytoplankton  bio¬ 
mass  in  the  pelagic  tropical  North  Atlantic,  and  0.2  to 
0.8%  of  daily  carbon  flux  from  the  global  ocean. 
Moreover,  different  species  of  seagrass  may  play 
different  ecological  roles  based  on  whether  the  major¬ 
ity  of  the  biomass  accumulates  locally  within  the  bed 
or  is  transported  out  of  the  ecosystem  (Zieman  et  al. 
1979,  Burdige  &  Zimmerman  2002,  Burdige  et  al. 
2008).  Hence  improved  methods  for  quantifying  large- 
scale  benthic  ecology  are  warranted. 

Quantifying  productivity  of  benthic  ecosystems  is 
challenging  due  to  the  high  degree  of  fine-scale  vari¬ 
ability  in  the  seafloor  features  in  many  environments 
(cm  to  m  scale)  and  the  limited  perspective  available 
from  traditional  in  situ  sampling  methods.  Remote 
sensing  of  benthic  environments  from  space  offers  the 
potential  for  filling  in  this  gap  (reviewed  in  Dekker  et 
al.  2006,  Phinn  et  al.  2008);  however,  such  methods 
require  considerable  in  situ  measurements  to  develop 
approaches  to  differentiate  among  types  of  benthic 
producers  (seagrass,  algal-covered  sand,  sediment) 
and  to  differentiate  the  optical  signatures  of  the  water 
column  and  the  benthos. 

The  GBB  is  an  ideal  experimental  platform  to  develop 
and  test  remote  sensing  approaches  for  estimating 
NPP  of  benthic  habitats.  The  water  is  consistently  shal¬ 
low  across  the  bank  (3  to  10  m)  and  large  regions  of 
'green'  seafloor  are  observable  from  space  (Fig.  1). 
Moreover,  the  water  overlying  much  of  the  GBB  is  rel¬ 
atively  oligotrophic,  without  significant  water  column 
phytoplankton  biomass.  Water  column  optical  proper¬ 
ties  are  primarily  dominated  colored  dissolved  organic 
matter  produced  by  benthic  organisms,  such  as  sea¬ 
grass  (Boss  &  Zaneveld  2003,  Otis  et  al.  2004).  How¬ 
ever,  on  the  western  portion  of  the  GBB  near  Andros 
Island,  fine  lime  mud  and  pellet  mud  sediments  can  be 
periodically  suspended  into  the  water  column,  produc¬ 


ing  ‘whitings’  that  are  visible  from  space  (Robbins  et 
al.  1997,  Dierssen  et  al.  2009a).  These  turbid,  wind-dri¬ 
ven  whiting  events  can  obscure  the  benthos  from  view, 
but  cover  only  a  relatively  small  area  per  episode. 
Hence,  the  expansive  GBB  is,  for  the  most  part,  opti¬ 
cally  shallow,  in  that  reflectance  of  the  seafloor  con¬ 
tributes  to  the  reflected  light  measured  by  a  satellite  or 
aircraft. 

As  with  phytoplankton  in  the  open  ocean,  variations 
in  benthic  environments  change  the  ocean  color 
expressed  at  the  sea  surface,  which  can  be  exploited  to 
estimate  the  distributions  of  the  respective  seafloor 
constituents.  Many  methods  have  been  developed  for 
conducting  remote  sensing  of  seagrasses  and  other 
benthic  producers  (reviewed  by  Kelly  1980,  Dekker  et 
al.  2006).  Such  approaches  vary  depending  on  the 
spectral  resolution  of  available  imagery  and  can  in¬ 
volve  look-up  tables  (Dierssen  et  al.  2003,  Louchard  et 
al.  2003,  Mobley  et  al.  2005),  hyperspectral  optimiza¬ 
tion  techniques  (Lee  et  al.  1999),  and  other  radiative 
transfer-based  approaches  (Dekker  et  al.  2005).  The 
amount  and  spectral  shape  of  reflected  light  is  influ¬ 
enced  by  water  depth,  water  column  optical  proper¬ 
ties,  and  reflectance  signatures  of  the  underlying 
benthos  (Armstrong  1993,  Maritorena  et  al.  1994). 
However,  methods  depend  on  the  spectral  resolution 
available  from  a  sensor  (Phinn  et  al.  2008),  in  addition 
to  the  spatial  resolution,  and  the  ability  to  retrieve  cal¬ 
ibrated  and  atmospherically  corrected  imagery  over 
the  optically  shallow  region. 

In  the  present  study,  field  data  were  used  to  develop 
an  ocean  color  remote  sensing  model  to  characterize 
the  benthic  habitats  of  the  GBB.  Using  a  radiative 
transfer-based  approach,  the  distribution  of  broad  cat¬ 
egories  of  primary  producer  habitats  (seagrass,  algae, 
and  sand)  were  estimated  from  variations  in  remotely 
sensed  ocean  color,  from  hues  of  light  blue  (primarily 
sand)  to  green  (vegetated),  across  the  GBB  and  com¬ 
bined  with  empirical  relationships  between  primary 
producer  abundance  and  benthic  productivity  to  esti¬ 
mate  the  contribution  of  these  highly  productive  shal¬ 
low  seas  to  the  global  carbon  cycle. 

MATERIALS  AND  METHODS 

The  present  study  was  conducted  on  the  GBB  in 
March  2004  aboard  the  RV  'Walton  Smith’  (Fig.  1). 
Twenty-five  stations  were  sampled  across  the  northern 
portions  of  GBB  near  Andros  and  along  the  Exumas 
Archipelago  (Fig.  1).  At  each  station,  divers  were 
deployed  in  the  water  to  survey  the  site,  quantify  sea¬ 
grass  densities  if  present,  and  collect  sediment  cores. 
Water  column  optical  and  physical  measurements  and 
hyperspectral  reflectance  measurements  were  con- 
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Fig.  1.  SeaWiFS  true  color  image  (6  March  2004)  of  the  tropical  West  Atlantic 
reveals  extensive  areas  of  green  and  blue-colored  optically  shallow  water 
across  the  entire  Great  Bahama  Banks.  Dominant  benthic  type  at  each  field 
station  is  identified  in  the  key 


ducted  in  situ.  Sediment,  seagrass,  and  macroalgae 
samples  were  collected  by  divers  and  analyzed  further 
in  the  laboratory.  Reflectance  measurements  were  con¬ 
ducted  on  the  sediment  cores  collected  at  all  stations. 
Oxygenic  photosynthesis  and  respiration  was  mea¬ 
sured  at  a  subset  of  the  stations  representing  the  differ¬ 
ent  habitats:  turtle  grass  leaves  of  all  ages  present 
within  a  shoot  were  analyzed  with  epiphytes  from 
Stns  WS12,  WS13,  WS16,  and  WS20;  the  macroalgae 
Colpomenia  sp.  was  analyzed  from  Stn  WS7  where  it 
was  found  in  organized  benthic  windrows  caused  by 
Langmuir  circulation  (Dierssen  et  al.  2009b);  grape- 
stone  sediments  with  extensive  microalgal  cover  were 
analyzed  from  Stn  WS6.  Seagrass  productivity  was  also 
estimated  from  in  situ  leaf  growth  rates  and  compared 
to  field  measurements  conducted  in  2002  and  2003 
near  Lee  Stocking  Island,  Bahamas  (Dierssen  et  al. 
2003). 


Optical  measurements.  A  variety  of  opti¬ 
cal  measurements  were  conducted  to  char¬ 
acterize  how  the  water  column  and  sea¬ 
floor  absorb,  scatter,  and  reflect  light.  The 
reader  is  referred  to  Mobley  (1994)  for  a 
more  detailed  characterization  of  optical 
terminology.  An  ac-9  Plus  package  (WET 
Labs)  was  deployed  at  each  station  to  mea¬ 
sure  light  absorption  and  attenuation  of  the 
water  column  at  a  fixed  depth  of  2  to  3  m 
below  the  sea  surface.  Salinity  and  temper¬ 
ature  were  measured  using  an  SBE  19 
(Sea-Bird  Electronics)  integrated  into  the 
package.  The  particulate  and  dissolved 
absorption  and  attenuation  coefficients,  apg 
and  cpg,  respectively,  were  corrected  for 
temperature,  salinity,  and  scattering  effects 
(Dierssen  et  al.  2009a).  Drift  offsets  were 
obtained  from  a  post-cruise  clean  water 
calibration.  Total  particulate  scattering  (bp) 
was  calculated  as  the  difference  between 
cpg  and  apg,  assuming  negligible  scattering 
from  dissolved  material.  The  particulate 
backscattering  coefficient  (fobp)  refers  to  all 
the  photons  that  have  been  redirected  in 
the  backward  direction  due  to  scattering 
from  particles  in  the  water  and,  to  a  first 
order,  is  positively  correlated  to  the  total 
concentration  of  particles  in  the  water  (re¬ 
viewed  in  Stramski  et  al.  2004).  The  ECO- 
VSF  (WET  Labs)  provided  backscattering 
measurements  at  3  angles  (100,  125,  and 
150°)  and  3  wavelengths  (470,  532,  and 
650  nm).  Data  were  corrected  for  light 
attenuation  effects  using  concurrent  ac-9 
data  and  extrapolated  from  90  to  180°  using 
a  third-order  polynomial.  The  backscatter¬ 
ing  ratio  (bbp)  was  calculated  as  the  ratio  of  btjP  to  bp. 

Remote  sensing  reflectance  (Rrs)  represents  the  spec¬ 
tral  distribution  of  light  emerging  from  the  sea  surface 
normalized  to  that  incident  upon  the  sea  surface  and 
was  measured  at  each  station  using  a  Field  Spec  Pro™ 
VNIR-NIR1  portable  spectrometer  system  (Analytical 
Spectral  Devices).  A  sequence  of  measurements  was 
made  with  an  8°  foreoptic  focused  at  a  40  to  45°  angle 
sequentially  on  a  gray  plaque,  sea  surface,  and  sky. 
The  parameter  RIS  was  calculated  as  the  water-leaving 
radiance  normalized  to  downwelling  irradiance  reach¬ 
ing  the  sea  surface  following  Mobley  et  al.  (1999).  Val¬ 
ues  of  RIS  exclude  photons  which  have  reflected  off  the 
sea  surface  and  do  not  penetrate  the  water  surface  (i.e. 
sun  glint).  Residual  reflected  sky  radiance  was  cor¬ 
rected  using  reflectance  differences  between  715  and 
735  nm,  an  approach  specifically  developed  for  coastal 
waters  (Gould  et  al.  2001). 
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Seagrass  morphology,  density,  and  growth.  Shoot 
density  of  the  seagrasses  Thalassia  testudinum  (turtle 
grass),  Syringodium  filiforme  (manatee  grass),  and 
Halodule  wrightii  (shoal  grass)  were  estimated  at  each 
station  using  15  to  20  randomly  located  0.04  m2 
quadrats  within  a  20  m  radius  around  a  central  buoy, 
spanning  roughly  1000  m2  at  each  station.  One  shoot  of 
each  species  (if  present)  was  collected  randomly  from 
each  quadrat  and  the  length  (nearest  mm)  and  width 
(nearest  0.01  mm)  of  each  leaf  was  measured  in  the 
laboratory.  Leaf  area  index  (LAI)  was  estimated  from 
leaf  morphometries  and  shoot  density  counts  as  the 
average  one-sided  leaf  area  per  surface  area  (m2)  of 
substrate  (Zimmerman  2003). 

Reflectance  of  individual  turtle  grass  leaves  were 
measured  in  the  laboratory  using  a  Shimadzu  UV2101- 
PC  scanning  spectrophotometer  fitted  with  an  inte¬ 
grating  sphere.  Sediment  cores  were  collected  by 
divers  and  the  surface  reflectance  was  estimated  in  the 
laboratory  using  the  Field  Spec-Pro  with  the  bare  fiber 
optic  probe  positioned  over  the  illuminated  sample  and 
referenced  against  a  white  spectralon  plaque  (99.9% 
reflectance)  positioned  similarly.  For  each  core,  5  dif¬ 
ferent  spectra  were  recorded  and  averaged.  Because 
of  the  way  light  refracts  as  it  transfers  through  an 
air-water  interface,  these  measurements  approximate 
the  spectral  shape  of  reflectance,  but  the  absolute 
spectral  magnitude  may  vary  somewhat  from  the 
reflectance  of  submerged  cores. 

Leaf  growth  rates  were  measured  at  2  stations  popu¬ 
lated  by  turtle  grass.  One  shoot  from  each  quadrat  was 
marked  with  a  20  gauge  hypodermic  needle  at  the  top 
of  the  leaf  sheath  (Zieman  &  Wetzel  1980)  and  tagged 
with  a  length  of  pink  flagging  tape  for  subsequent 
identification.  The  shoots  were  collected  2  wk  after 
marking  and  the  new  growth  of  each  leaf  was  mea¬ 
sured.  NPP  was  estimated  for  each  site  as  the  product 
of  LAI  and  leaf  growth  rate  (%  d_1)  and  converted  to  a 
mass-specific  rate  using  measured  ratios  for  shoot 
fresh  weight  (FW)  and  dry  weight  (DW),  which  ranged 
from  0.09  to  0.22  g  DW  g-1  FW.  Average  conversions  of 
0.05  g  FW  cm"2,  0.16  g  DW  g"1  FW,  and  0.35  gC  g"1  DW 
were  used  (Zimmerman  2003). 

Primary  producer  metabolism.  Oxygenic  photosyn¬ 
thesis  and  respiration  of  turtle  grass  leaves,  brown  al¬ 
gae  Colpomenia  sp.  thalli,  and  grapestone  sediments 
were  measured  polarographically  in  clear  glass  water- 
jacketed  incubation  chambers  (5  ml  volume).  For  the 
turtle  grass  leaves,  3  replicates  of  the  different  leaf  ages 
were  conducted  at  Stn  WS04,  2  replicates  at  Stns  WS12 
and  WS13,  and  1  measurement  for  each  leaf  age  at 
Stns  WS02,  WS03,  WS16,  and  WS20.  Three  replicates 
were  conducted  for  Colpomenia  sp.  and  only  1  mea¬ 
surement  for  the  grapestone  sediment.  Temperature  in¬ 
side  the  chambers  was  maintained  at  22°C  and  elec¬ 


trode  drift  was  measured  in  chambers  filled  with  air- 
saturated  seawater  prior  to  the  addition  of  each  sample. 
Respiration  of  each  sample  was  measured  in  the  dark 
using  air-saturated  seawater.  At  the  end  of  the  respira¬ 
tion  period,  02  content  of  the  incubation  medium  was 
further  reduced  to  about  50%  of  air  saturation  by  bub¬ 
bling  with  a  gas  mixture  of  C02  (350  ppm)  and  N2  (bal¬ 
ance)  before  commencing  photosynthesis  measure¬ 
ments  The  light-saturated  rate  of  photosynthesis  (Pm) 
was  measured  by  exposing  each  sample  to  -600  pmol 
quanta  m-2  s_1  of  white  light  provided  by  a  fiber  optic 
light  source  using  a  variable  intensity  incandescent 
lamp  (Dolan-Jenner  Industries,  Fiber-Lite  Series  180). 
Incubation  irradiance  was  measured  using  a  factory- 
calibrated  quantum  irradiance  meter  (Biospherical  In¬ 
struments,  QSL-100).  Sample  weights  were  determined 
after  drying  at  65°C  for  -10  d.  NPP  was  estimated  as¬ 
suming  photosynthesis  was  light  saturated  for  12  h  d_1 
and  the  photosynthetic  quotient  (02:C02)  was  1.2, 
following  the  procedure  outlined  in  Dierssen  et  al. 
(2009b).  Seagrass  leaves  were  assumed  to  comprise 
0.05  g  FW  cm-2  of  leaf  area,  which  was  derived  from 
turtle  grass  samples  collected  at  Lee  Stocking  Island 
(Zimmerman  2003).  The  drifting  macroalga  Colpome¬ 
nia  sp.  was  determined  to  be  3  to  10  cm  in  thickness  dis¬ 
tributed  in  ordered  windrows  (Dierssen  et  al.  2009b) 
and  assumed  to  be  2  cm  in  thickness  averaged  over 
non-windrow  area  without  windrows.  The  macroalgae 
was  negatively  buoyant,  but  the  density  was  conserva¬ 
tively  assumed  to  be  that  of  seawater  (1025  kg  mr3). 

Image  processing.  A  full  resolution  (1  km)  Sea-view- 
ing  Wide  Field-of-View  (SeaWiFS)  Level  1A  Local 
Area  Coverage  image  from  6  March  2004  provided  a 
relatively  cloud-free  view  of  the  GBB.  Standard  atmos¬ 
pheric  correction  routines  for  SeaWiFS,  which  were 
developed  for  optically-deep  open  ocean  water,  cannot 
be  applied  accurately  to  the  optically  shallow  GBB. 
Consequently,  RIS  was  derived  at  each  pixel  in  the 
image  by  applying  an  iterative  scheme  to  first  process 
the  image  using  the  SeaWiFS  multiple-scattering 
atmospheric  correction  algorithm  (Gordon  1997).  The 
atmospheric  correction  scheme  from  neighboring  deep 
water  was  identified  and  the  entire  image  was  then 
reprocessed  with  SeaWiFS  Data  Analysis  System  using 
a  fixed  aerosol  model  for  the  region.  This  type  of 
atmospheric  correction  is  only  an  approximation  of  the 
atmospheric  properties  over  shallow  water.  Errors  in 
aerosols,  whitecaps,  and  sun  glint  estimated  over  deep 
water  can  lead  to  consistent  biases  in  retrieved 
reflectances  over  shallow  water.  A  linear  correction 
was  derived  for  each  wavelength  where  the  offset 
represents  errors  due  to  path  radiance  and  the  slope 
represents  errors  in  the  transmission  of  water-leaving 
radiance  through  the  atmosphere  (Gordon  1997). 
Linear  regression  was  conducted  with  RIS  derived  from 
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the  SeaWiFS  sensor  and  in  situ  measurements 
from  a  subset  of  6  stations  considered  to  be  consis¬ 
tent  at  km  scales  from  our  field  observations.  Some 
stations,  for  example,  had  patchy  distributions  of 
seagrass  and  were  excluded  from  this  analysis. 
The  image  was  then  georeferenced  and  warped 
based  on  a  complete  geometry  model  of  the  earth 
and  satellite  orbit.  Estimates  of  bottom  reflectance, 
Rb,  were  modeled  from  the  RIS  (in  situ  values  used 
for  model  development  or  satellite-derived  values 
for  model  application)  following  the  look-up  table 
approach  outlined  in  Dierssen  et  al.  (2003),  as  dis¬ 
cussed  further  in  the  'Results'. 


RESULTS 

Benthic  characterization 

We  encountered  a  variety  of  benthic  habitats 
across  the  study  area,  including  those  predomi¬ 
nantly  comprised  of  seagrass,  benthic  macroalgae, 
and  sediment  microalgae.  Seagrass  meadows 
(Fig.  1)  were  concentrated  along  the  outer  perimeter 
of  the  bank  in  areas  experiencing  regular  tidal  ex¬ 
change  with  Atlantic  waters.  Stations  to  the  west  of 
Andros  Island  were  comprised  of  fine  pellet  muds, 
but  also  contained  sparse  distributions  of  seagrass. 
Along  the  Exumas,  4  different  types  of  benthic  com¬ 
munities  were  found:  seagrass,  brown  macroalgae, 
red  macroalgae,  and  microalgae  associated  with 
grapestone  sediment.  Bahamian  grapestone  is  a 
coarse  sediment  produced  by  microbial  colonization 
and  cementation  of  recrystallized  ooids  that  gives 
the  aggregated  grains  the  appearance  of  a  bunch  of 
grapes  (Winland  &  Matthews  1974).  The  algal  con¬ 
tent  of  grapestone  biofilms  can  be  highly  variable, 
with  chlorophyll  concentrations  ranging  from  2  to 
10  pg  g_I  (Stephens  et  al.  2003).  At  stations  sampled 
along  the  Exumas,  the  grapestone  color  ranged  from 
light  to  dark  green,  indicative  of  variability  in  micro- 
algal  density.  The  phaeophyte  Colpomenia  sp.  was 
observed  at  the  westernmost  station  sampled  on  the 
Exumas  in  a  series  of  windrows  ordered  along  the 
seafloor  by  wind-driven  helical  (Langmuir)  circula¬ 
tion  (Monahan  1993,  Dierssen  et  al.  2009b). 

LAI  provided  a  quantitative  measure  of  the  sea¬ 
grass  abundance  that  incorporated  both  the  shoot 
density  and  leaf  morphometric  data  (height  and 
width).  The  highest  LAI  on  the  GBB  (3.7  ±  0.4)  was 
measured  at  a  station  near  Bimini  Island  (WS01, 
5.2  m  depth)  containing  a  mixture  of  the  cylindri- 
cally  bladed  manatee  grass  and  flat-bladed  turtle 
grass  (Table  1).  Other  dense  seagrass  meadows 
comprised  of  both  species  (LAI  1.5  to  2.2)  were 
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found  further  north  in  this  area  at  depths  up  to  10  m 
(Fig.  1).  Along  the  Exumas,  meadows  of  turtle  grass 
occurred  primarily  at  3  to  5  m  depth  in  water  bordering 
the  channels  near  the  islands  with  LAI  ranging  up  to  2. 
At  2  stations  close  to  the  Exuma  Islands,  shoal  grass 
was  found  interspersed  with  turtle  grass.  Sparse  distri¬ 
butions  of  turtle  grass  were  found  directly  west  of 
Andros  Island  growing  in  the  mud  sediment  (LAI  <  0.2). 
Isolated  patches  of  dense  turtle  grass  occurred  sporad¬ 
ically  along  Andros  (WS22)  surrounded  by  extensive 
flats  of  white  sand,  sponges,  and  isolated  coral  heads. 

Total  LAI  for  all  species  of  seagrass  could  be  pre¬ 
dicted  by  total  shoot  density  (shoots  m-2),  which  ex¬ 
plained  over  80  %  of  the  variability  in  the  data  (Fig.  2A). 
Although  the  regression  was  heavily  influenced  by  one 
station  with  very  high  LAI,  the  general  relationship 
between  shoot  density  and  LAI  was  strongly  linear, 


Fig.  2.  (A)  Relationship  between  seagrass  shoot  density  and 
estimated  leaf  area  index  (LAI)  ranging  from  0  to  3.5.  (B)  Sea¬ 
grass  LAI  and  above-ground  net  primary  production  (NPP) 
determined  in  2002-2003  near  Lee  Stocking  Island,  Bahamas, 
from  in  situ  leaf  growth  rates  (o),  in  2004  at  northern  Bimini 
stations  from  in  situ  growth  rates  (A),  and  in  2004  from  labora¬ 
tory  metabolism  measurements  (A).  Dashed  lines  represent 
95  %  confidence  intervals  for  the  regression 


despite  the  fact  that  different  stations  contained  differ¬ 
ent  mixtures  of  seagrass  species.  The  residual  variabil¬ 
ity  among  stations  was  due  to  differences  in  shoot  mor¬ 
phometry,  including  leaf  length  and  width.  For  both 
seagrass  species,  canopy  heights  were  significantly 
taller  at  the  Andros  stations  compared  to  the  Exumas. 
Shoot  leaf  area  of  turtle  grass,  calculated  as  the  sum  of 
the  one-sided  area  of  all  leaves  on  each  shoot,  averaged 
29  cm2  shoot-1  for  Andros  stations  and  only  15  cm2 
shoot-1  for  the  Exumas  populations.  Shoot  leaf  area  of 
manatee  grass  was  11  cm2  shoot-1  at  Andros  stations 
and  only  3  cm2  shoot-1  at  the  Exumas.  No  trend  was 
evident  between  the  canopy  height  and  depth  of 
occurrence  (r2  =  0.17,  p  >  0.05,  data  not  shown). 

The  ratio  of  gross  photosynthesis  to  respiration  for 
turtle  grass  averaged  5.19  ±  2.74  (Table  2),  resulting  in 
NPP  ranging  from  0.2  ±  0.1  gC  m-2  d-1  for  relatively 
sparse  seagrass  with  an  LAI  of  0.5  up  to  0.8  ±  0.3  gC 
m-2  d-1  for  dense  seagrass  with  an  LAI  of  2.5.  NPP  from 
in  situ  leaf  growth  rates  from  2  Andros  stations  (WS03 
and  WS04)  measured  over  a  2  wk  period  were  similar 
in  magnitude  (Fig.  2B,  solid  triangles)  and  followed  the 
relationship:  NPP  =  0.4  x  LAI. 

In  comparison,  NPP  estimated  in  May-June  2002 
and  2003  near  Lee  Stocking  Island  on  the  Exumas 
(Dierssen  et  al.  2003)  was  nearly  2-fold  higher  (slope  of 
0.78)  than  2004  rates  (Fig.  2B,  circles).  This  higher  rate 
could  be  due  to  differences  between  spring  and  sum¬ 
mer  sampling  and  greater  exchange  of  nutrients  on  the 
Exumas  compared  to  Andros. 

Rates  of  NPP  for  benthic  algae  were  quite  variable  de¬ 
pending  on  the  taxon.  The  highest  metabolism  and  cor¬ 
responding  NPP  was  measured  for  the  phaeophyte  Col- 
pomenia  sp.  found  growing  rapidly  on  the  mid-portion 
of  the  Exumas.  Dense  accumulations  of  this  brown 
macroalga  were  found  as  unattached  free-drifting  mats 
tumbling  along  the  bottom  with  currents  believed  to  be 
generated  by  Langmuir  circulation  (Dierssen  et  al. 
2009b).  From  photosynthesis  estimates  on  individual 
thalli,  we  estimated  NPP  to  be  18  ±  4.4  gC  m-2  d-1.  How¬ 
ever,  this  estimate  does  not  incorporate  self-shading  and 
spectral  changes  in  light  availability  within  the  drifting 
mats  (Vergara  et  al.  1997,  Brush  &  Nixon  2003),  which 
may  serve  to  diminish  NPP  somewhat. 

In  contrast,  grapestone  sediments  containing  thin 
algal  biofilms  were  barely  net  autotrophic,  with  rates 
of  0.005  gC  m-2  d-1.  NPP  of  red  algae  and  grapestone 
with  dense  biofilm  was  not  measured  directly.  Benthic 
microalgae  within  sediment  from  the  seagrass  mead¬ 
ows  may  also  contribute  to  NPP  (MacIntyre  et  al.  1996) 
but  was  not  measured  directly  in  the  present  study. 
Around  Lee  Stocking  Island  Bahamas,  sediment 
chlorophyll  a  (chi  a)  concentrations  were  3-fold  higher 
in  dense  turtle  grass  beds  (3  to  4  pg  g-1  sediment)  com¬ 
pared  to  adjacent  sediment  (Stephens  et  al.  2003). 
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Table  2.  Thalassia  testudinum,  Colpomenia  sp.  and  grapestone  sediment.  Net  primary  productivity  (NPP)  of  benthic  producers 
measured  across  the  Great  Bahama  Bank.  Leaf  area  index  (LAI)  is  the  area  (m2)  of  1 -sided  seagrass  leaf  per  m2  seafloor. 

Pg:R:  ratio  of  gross  production  to  respiration;  DW:  dry  weight 


Producer 

Gross  oxygen  evolution,  Pg 
(pmol  02  min-1  g  DW-1) 

Pg:R 

Areal  average 
(g  DW  nr2) 

NPP 

(gC  nr2  d-1) 

Seagrass  Thalassia  testudinum 

0.78  ±  0.28  (n  =  23) 

5.19  ±  2.74 

LAI  0.5 

50 

0.17  ±  0.11 

LAI  1.0 

100 

0.35  ±  0.13 

LAI  1.5 

150 

0.52  ±  0.19 

LAI  2.0 

200 

0.69  ±  0.25 

Macroalgae  Colpomenia  sp. 

2.20  ±  0.24  (n  =  3) 

4.71  ±  0.59 

2050 

18  ±4.4“ 

Grapestone  algae 

0.056  (n  =  1) 

2.12 

205 

0.005 

“Assuming  no  self-shading  or  changes  in  spectral  composition  of  light  within  the  2  cm  thick  drifting  macroalgal  mats 

Optical  properties 

In  order  to  interpret  a  remote  sensing  signal  derived 
from  the  benthic  primary  producers,  water  column  opti¬ 
cal  properties  were  measured  in  the  surface  waters  at 
all  benthic  stations,  as  well  as  in  surface  waters  of  the 
Tongue  of  the  Ocean  (>1000  m).  Consistent  with  past 
studies  and  visual  observations  by  the  divers,  the  absorp¬ 
tion  (apg)  and  attenuation  (cpg)  coefficients  were  gener¬ 
ally  low  across  the  entire  region,  indicating  high  water 
clarity  (Fig.  3A,B).  The  spectral  shape  of  the  absorption 
coefficient  was  consistent  with  the  presence  of  colored 
dissolved  organic  matter  and  few  phytoplankton  sus¬ 


pended  in  the  water  column  (Boss  &  Zaneveld  2003). 
Water  from  the  Tongue  had  considerably  lower  absorp¬ 
tion  and  attenuation  coefficients  compared  to  the  bank. 
Andros  also  had  higher  cpg  than  either  the  Exumas  or  the 
Tongue,  likely  due  to  the  presence  of  non-absorbing 
'white'  calcite  fine  sediment  particles  suspended  in  the 
water  column  (Shinn  et  al.  1989,  Dierssen  et  al.  2009a). 
We  assumed  that  apg  was  primarily  due  to  dissolved  ab¬ 
sorption  (ag)  and  subtracted  this  from  the  cpg  values  to 
obtain  an  estimate  of  light  attenuation  of  particulates 
alone,  cp.  The  hyperbolic  wavelength  spectra  of  cp  for 
Andros,  Exumas,  and  the  Tongue,  respectively,  were 
0.95  ±  0.21,  0.98  ±  0.22,  and  0.8  (1  measurement  only). 


Fig.  3.  Mean  inherent 
optical  properties  (±1 
SD)  measured  at  sta¬ 
tions  from  Andros,  Exu¬ 
mas  and  the  Tongue  of 
the  Ocean:  (A)  absorp¬ 
tion  by  particulate  and 
dissolved  material  (apg); 
(B)  attenuation  by  par¬ 
ticulate  and  dissolved 
material  (cpg);  (C)  par¬ 
ticulate  backscattering 
(bbp);  and  (D)  backscat¬ 
tering  ratio  (bbp) 


Wavelength  (nm) 
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The  particulate  backscattering  data  revealed  a  simi¬ 
lar  pattern.  Levels  of  p  from  stations  west  of  Andros 
were  twice  as  high  as  stations  near  the  Exumas 
(Fig.  3C),  indicative  of  more  suspended  sediment  on 
the  western  arm  of  the  GBB.  This  difference  was  likely 
due  to  the  larger  sediment  grains  found  along  the  Exu¬ 
mas  which  were  less  prone  to  resuspension  (Hu  2007, 
Dierssen  et  al.  2009a).  Water  from  the  Tongue  was  so 
devoid  of  particles  that  measurements  of  were 
below  the  limit  of  detection  for  the  instrument  (data  not 
shown).  The  ratio  of  backscattering  to  total  scattering 
(hbp)  was  high  and  consistent  between  the  Andros  and 
Exumas  regions  (Fig.  3D).  Measurements  averaged 
-0.028  across  all  stations  and  no  significant  wave¬ 
length  dependence  was  observed.  Such  high  values 
are  consistent  with  the  presence  of  minerogenic  parti¬ 


cles  (i.e.  suspended  sediment)  rather  than  biogenic 
particles  (i.e.  phytoplankton).  Levels  of  fobp  reached  as 
high  as  0.05  during  an  intense  resuspension  events  or 
'whiting'  which  occurred  at  WS19  west  of  Andros 
Island  (Dierssen  et  al.  2009a).  The  estimated  bulk 
index  of  refraction  calculated  according  to  Twardow- 
ski  et  al.  (2001)  on  the  GBB  was  -1.16,  consistent  with 
the  presence  of  inorganic  minerogenic  sediment  parti¬ 
cles  and  few  organic  phytoplankton. 

Sea  spectral  reflectance  used  for  remote  sensing  stud¬ 
ies  {RIS,  i.e.  water-leaving  radiance  normalized  to  inci¬ 
dent  downwelling  irradiance)  provided  an  indicator  of 
the  water  color  across  the  visible  spectrum  at  each  sta¬ 
tion.  Levels  of  RIS  from  selected  locations  at  5  m  depth 
revealed  quite  different  spectral  magnitudes  and  shapes 
depending  on  the  seafloor  composition  (Fig.  4A).  Grape- 


Organic  sand,  5.2  m 

— •- 

-  Andros  mud,  5.4  m 

Grapestone  light,  5.6  m 

Sparse  seagrass,  5.5  m 

-  Red  macroalgae,  5.6  m 

Brown  macroalgae,  5.3  m 

-  Dense  seagrass,  5.3  m 

- -x— 

-  Grapestone  dark,  5.7  m 

Fig.  4.  (A)  Spectral  remote  sensing  reflectance 
(Rrs)  measured  at  8  stations  with  -5  m  water 
depth  showing  the  decrease  in  magnitude 
from  bright  sand  to  dark  green  grapestone 
sediment  (sr  =  steradian).  (B)  The  same  spec¬ 
trum  from  (A)  normalized  to  Rrs(555)  showing 
the  spectral-shift  from  blue-  to  green-peaked 
with  more  productive  benthos.  (C)  Derivative 
spectrum  of  RIS  showing  wavebands  with 
differences  in  spectral  slope.  (D)  Remote  sens¬ 
ing  reflectance  over  dense  seagrass  beds 
ranging  in  depth  from  5  to  10  m.  (E)  Labora¬ 
tory  measurements  of  benthic  reflectance 
(Rb)  for  different  types  of  sediment,  seagrass 
leaves,  and  macroalgal  blades 
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stone  sediment  with  a  green  biofilm  had  the  lowest  RIS 
and  peaked  in  the  green  at  -570  nm.  Stations  com¬ 
prised  of  red  and  brown  macroalgae  had  slightly 
higher  RIS  but  similar  spectral  shapes.  Dense  seagrass 
beds  had  a  similar  magnitude  in  the  green  (500  to 
570  nm),  but  appeared  to  reflect  more  light  in  the  blue 
portion  of  the  spectrum  (400  to  500  nm)  than  the  other 
benthic  producers.  Remotely  sensed  reflectance  from 
organic  sand  and  Andros  mud  were  both  -5  times 
greater  than  RIS  from  the  productive  benthic  habitats, 
including  grapestone  sediment  with  microalgal  bio¬ 
film.  Moreover,  the  spectra  peaked  at  500  nm  com¬ 
pared  to  the  550  to  570  nm  peak  observed  for  the  other 
spectra. 

When  the  spectra  were  normalized  to  RIS  at  555  nm, 
shape  differences  between  the  various  benthic  types 
were  readily  observed  (Fig.  4B).  At  5  m  depth,  the 
reflectance  peaks  of  high  productivity  stations  were 
green-shifted  relative  to  low  productivity  stations. 
Fourth  derivatives  of  the  spectra  (Fig.  4C)  revealed 
spectral  regions  that  could  potentially  be  used  to  dif¬ 
ferentiate  different  benthic  types  with  high  signal-to- 
noise  imaging  spectroscopy.  Subtleties  in  spectral 
shape  were  most  apparent  in  the  green  region  of  the 
spectrum  with  differences  observed  at  505  to  510,  515, 
525,  and  550  nm. 

Highly  productive  dense  seagrass  meadows  were 
found  growing  at  depths  >10  m  at  the  north  end  of 
Andros  (WS03  and  WS04)  and  the  reflectance  signa¬ 
ture  was  dominated  by  blue  light  rather  than  green  as 
the  influence  of  water  column  on  RIS  increased  with 
depth  (Fig.  4D).  From  space,  therefore,  the  deep  sea¬ 
grass  beds  appeared  more  blue  than  green.  In  addi¬ 
tion,  reflectance  in  the  red  portion  of  the  spectrum  (600 
to  700  nm)  was  also  2-  to  3-fold  greater  at  the  shallower 
stations,  because  more  of  the  red  light  reflected  from 
seafloor  reaches  the  sea  surface  (Dierssen  et  al.  2003). 

The  reflectance  of  the  benthos  (RB)  associated  with 
each  station  was  a  combination  of  differently  colored 
sediments  and  benthic  constituents.  Reflectance  from 
seagrass  leaves  and  macroalgae  blades  measured  in 
the  laboratory  (Fig.  4E)  had  low  reflectance  across 
most  of  the  visible  spectrum.  Seagrass  leaves  were 
more  reflective  in  the  blue  and  green  compared  to 
brown  macroalgae,  which  peaked  further  in  the  red  at 
600  nm.  However,  the  seagrass  leaf  spectrum  does  not 
represent  the  reflectance  from  the  seagrass  canopy, 
which  is  a  spectral  mixture  of  darker  seagrass  leaves 
and  brighter  sediment  that  varies  depending  on  the 
density  of  the  canopy  and  the  brightness  of  the  under¬ 
lying  sand  (Zimmerman  2003).  Sediment  RIS  from  dif¬ 
ferent  parts  of  the  GBB  had  a  large  range  in  bottom 
reflectance.  Measured  RB  from  grapestone  covered 
with  dense  algal  biofilm  was  low  in  magnitude  indicat¬ 
ing  considerable  pigment  concentrations.  In  contrast, 


sediment  from  the  Exumas  was  twice  as  reflective 
across  the  visible  spectrum.  The  fine  mud  sediment 
from  the  stations  around  Andros  was  even  more  reflec¬ 
tive,  peaking  at  around  0.4  between  550  and  650  nm. 
Such  large  differences  challenge  our  ability  to  estimate 
benthic  properties  from  multispectral  imagery,  if  bathy¬ 
metry  is  unquantified. 

Quantifying  basin-wide  productivity 

Our  basic  approach  for  modeling  NPP  involved  esti¬ 
mating  Rb  from  the  imagery  and  then  quantifying  pro¬ 
ductivity.  To  account  for  absorption  and  scattering 
within  the  water  column,  RB  was  estimated  from  RIS 
using  a  depth-specific  look-up  table  derived  from 
radiative  transfer  modeling  (Hydrolight,  Sequoia  Sci¬ 
entific;  Mobley  1994)  with  the  measured  optical  prop¬ 
erties  of  the  water  column  (Dierssen  et  al.  2003).  Aver¬ 
age  water  column  optical  properties  measured  during 
the  field  campaign  (Fig.  3,  Exumas)  were  used  in  the 
model  and  assumed  to  be  generally  representative  of 
the  GBB.  Fig.  5A  presents  the  blue/green  sea  surface 
reflectance  ratio  (Rrs[490]:Rrs[555J),  used  to  estimate 
chlorophyll  in  the  optically  deep  majority  of  the  world's 
oceans,  for  stations  comprised  of  different  benthic  con¬ 
stituents  and  at  different  depths.  Many  of  the  benthic 
stations  had  lower  reflective  ratios  indicative  of  chloro¬ 
phyll-containing  compounds  that  absorb  blue  light. 
However,  reflectance  ratios  of  the  deeper  seagrass  sta¬ 
tions  (9  to  10  m  depth)  mapped  into  the  blue  (less 
productive)  region  because  of  the  impact  of  the  over- 
lying  water  column.  When  the  spectral  influence  of  the 
overlying  water  column  was  removed  by  applying 
the  look-up  table  approach,  seagrass,  macroalgae,  and 
grapestone  sediment  with  dense  biofilm  all  grouped 
together  with  low  green  benthic  reflectance  RB(555), 
indicative  of  high  light  absorption,  while  less  produc¬ 
tive  stations  had  higher  modeled  RB(555)  (Fig.  5B). 

A  logarithmic  relationship  was  developed  from  the 
in  situ  data  to  approximate  NPP  over  2  orders  of  mag¬ 
nitude  from  modeled  RB(555)  (Fig.  5C,  R2  =  0.49,  p  < 
0.01,  n  =  26):  log(NPP)  =  -7.58RB(555)  +  0.65. 

The  relationship  between  benthic  reflectance  and 
NPP  separated  into  3  general  groups:  high,  medium, 
and  low  NPP.  High  NPP  was  associated  with  a  light 
absorbing  bottom  (i.e.  dark  reflectance);  low  to  non- 
detectable  NPP  was  associated  with  low  absorbing; 
and  bright  bottom  reflectances  (>3%)  were  character¬ 
istic  of  unvegetated  sands.  The  NPP  for  brown  macro¬ 
algae  Colpomenia  sp.  was  considerably  higher  than 
other  producers,  due  in  part  to  its  considerable  density 
on  the  seafloor  (Dierssen  et  al.  2009b).  NPP  for  red 
algae  was  assumed  to  be  similar  to  brown  algae,  but 
the  density  of  red  algae  was  estimated  by  the  divers  to 
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Fig.  5.  (A)  Ratio  of  blue  to  green  reflectance  was  higher  for  deeper  seagrass 
beds  compared  to  similar  density  beds  at  shallower  depths.  (B)  Modeled  green 
bottom  reflectance,  i?B(555),  from  the  in  situ  data  using  a  radiative  transfer- 
based  look-up  table.  (C)  An  inverse  relationship  was  found  between  the  mod¬ 
eled  Rb(555)  and  logarithmic  net  primary  productivity  (NPP)  differentiating 
regions  with  high,  medium,  and  low  NPP.  Dashed  lines  represent  95  %  confi¬ 
dence  intervals  for  the  regression 


be  at  least  1  order  of  magnitude  lower  than  Colpome- 
nia  sp.  For  the  stations  with  grapestone  sediment  cov¬ 
ered  with  dense  biofilm,  NPP  rates  were  estimated  to 
be  5  times  the  NPP  measured  for  the  thin  biofilm 
grapestone.  Past  studies  in  this  region  have  shown 
chi  a  concentrations  to  be  roughly  5  times  greater  in 
dense  biofilm  grapestone  compared  to  thin  biofilm 
grapestone  (10.2  versus  1.8  pg  g_1;  Stephens  et  al. 
2003).  The  present  study  was  limited  to  NPP  measure¬ 
ments  of  micro-  and  macroalgae  at  only  a  few  stations 
and  more  comprehensive  sampling  would  be  desir¬ 
able.  However,  we  feel  that  this  approach  represents 
the  general  relationship  between  the  greenness  of  the 
seafloor  and  the  amount  of  productivity,  and  captures 
the  large-scale  NPP  features.  The  majority  of  the  data 
points  (13  of  22)  fell  within  the  95%  confidence  inter¬ 
vals  and  within  a  factor  of  2  of  the  measurements. 

The  logarithmic  relationship  was  applied  to  SeaWiFS 
imagery  from  6  March  2004,  a  relatively  clear  sky 
satellite  overpass  from  1  wk  prior  to  the  start  of  the 
held  campaign.  Modeling  of  RB( 555)  required  accurate 
knowledge  of  the  bathymetry  in  order  to  access  the 


depth-dependent  look-up  table.  Sound¬ 
ings  across  the  bank  were  digitized 
(Maptech  2004)  and  gridded  to  the 
same  pixel  size  as  the  SeaWiFS  imagery 
(Fig.  6A).  Deeper  water  (>8  m)  occurred 
to  the  north  of  Andros  and  on  unsampled 
regions  on  the  southwest  of  the  GBB. 
Shallow  regions  occurred  to  the  west  of 
Andros  Island  and  along  much  of  the 
Exumas  (<5  m).  The  bathymetry  was 
corrected  for  tidal  height  to  match  bathy¬ 
metry  measured  at  the  sampling  stations. 
Accurate  application  of  the  NPP  relation¬ 
ship  also  necessitated  appropriate  atmo¬ 
spheric  correction  of  the  imagery.  As 
shown  in  Fig.  6B,  retrieved  estimates  of 
RIS  were  similar  in  both  magnitude  and 
spectral  shape  to  measurements  made 
at  stations  sampled  during  the  field 
campaign.  Values  of  SeaWiFS-derived 
RIS(5 55)  were  within  6%  of  the  measured 
spectrum  of  the  dark  target  and  18%  of 
the  bright  target. 

Satellite-derived  NPP  across  the  GBB 
(Fig.  6C)  was  low  in  regions  west  of  An¬ 
dros  Island  and  high  along  its  northern 
border.  Modeled  NPP  was  in  the  range  of 
held  measurements  from  stations  through¬ 
out  the  bank  and  generally  followed  a  1 : 1 
line  (r2  =  0.64,  p  <  0.01,  slope  =  1.2,  inter¬ 
cept  =  0.02;  Fig.  6D).  The  95  %  confidence 
intervals  fell  within  a  factor  of  1.7  to  2.8 
over  the  range  of  NPP  measured.  The  in 
situ  held  stations  were  sampled  at  a  spatial  scale  of 
-1000  m2  (assuming  divers  sampled  the  seafloor  ran¬ 
domly  within  a  20  m  radius  around  a  central  buoy)  com¬ 
pared  to  the  1  x  106  m2  area  of  an  image  pixel,  and  the 
validation  included  only  those  stations  with  measurable 
NPP  (i.e.  not  bare  sediment)  which  represented  large- 
scale  benthic  features  and  had  roughly  similar  bathym¬ 
etry  across  the  satellite  pixel  compared  to  the  measured 
depths.  The  median  NPP  across  the  GBB  was  0.38  gC 
m-2  d_1.  Total  NPP  aggregated  from  the  pixel  estimates 
across  the  95  000  km2  of  the  GBB  was  5.0  x  1010  gC  d-1. 
Assuming  these  rates  of  productivity  are  representative 
of  biomass  and  metabolic  activity  throughout  the  year, 
annual  NPP  would  be  on  the  order  of  1.8  x  1013  gC  yr_1. 

The  relative  area  covered  by  algae  (either  macro-  or 
microalgae)  versus  seagrass  (Fig.  7)  was  estimated  by 
evaluating  the  retrieved  bottom  reflectance  in  blue  and 
green  wavelengths  (i?B[440]  and  f?B[555]).  As  shown  in 
Fig.  4,  the  reflectance  magnitude  at  555  nm  was  simi¬ 
lar  for  dense  concentrations  of  all  3  types  of  benthic 
producers,  but  reflectance  in  the  blue  wavelengths 
(440  to  490  nm)  was  elevated  for  the  seagrass  mead- 
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Fig.  6.  (A)  Digitized  and  gridded  bathymetry  across  the  Great  Bahama  Bank  from  0  to  10  m.  (B)  Satellite-derived  remote  sensing 
reflectance  (RIS)  from  SeaWiFS  imagery  (O)  in  comparison  to  field  measurements  from  3  diverse  stations  (solid  lines),  sr  =  stera- 
dians  (C)  Modeled  net  primary  productivity  (NPP,  gC  nr2  yr-1)  derived  from  SeaWiFS  imagery  over  the  entire  Bahamas  Banks. 
(D)  Logarithmic  comparison  of  satellite-modeled  NPP  from  stations  throughout  the  banks  and  estimates  made  during  the  field 
campaign  (r2  =  0.64,  p  <  0.01).  Thin  dashed  lines  represent  95%  confidence  intervals  and  the  thick  dashed  line  is  the  1:1  line 


ows.  For  those  regions  with  NPP  >  0.3  gC  m-2  d"1,  we 
evaluated  the  retrieved  RB(AA0):RB(555)  ratio  derived 
from  the  imagery.  A  higher  ratio  was  presumed  to  be 
predominantly  seagrass  and  a  lower  ratio  was  either 
macro-  or  microalgae.  An  analysis  was  conducted  to 
determine  the  empirical  threshold  of  the  blue:green  RB 
ratio  which  best  discriminated  seagrass  from  the  algae. 
The  threshold  was  allowed  to  vary  in  0.01  increments 
from  0  to  2  until  the  highest  agreement  was  reached 
between  the  model  and  the  in  situ  identification  of  the 
benthos.  Seagrass  was  modeled  as  those  pixels  with  a 
blue: green  RB  ratio  >0.7  and  algae  below  this  empirical 
threshold  (Fig.  7A).  Sand  was  considered  to  be  pixels 
with  modeled  NPP  <0.3  gC  m-2  d_1. 

A  classification  accuracy  matrix  (Fig.  7B)  illustrates 
the  model  which  correctly  retrieved  the  highest  num¬ 
ber  of  stations.  The  rows  in  the  matrix  represent  the 
dominant  benthic  constituent  identified  at  the  field  sta¬ 
tions  and  the  columns  represent  the  benthic  class  iden¬ 
tified  from  the  satellite  imagery.  If  all  of  the  stations 
were  mapped  correctly,  then  numbers  will  only  appear 
in  the  diagonal  boxes  where  each  benthic  class  inter¬ 


sects;  in  all  cases,  the  sum  of  the  diagonals  is  used  to 
calculate  the  overall  accuracy.  Benthic  classes  at  17  of 
the  23  field  stations  (73%)  were  correctly  identified  by 
the  classification  technique.  The  off  diagonals  show 
errors  where  a  portion  of  the  map  is  omitted  from  the 
correct  category  (omission  error)  and  placed  in  an 
incorrect  category  (commission  error).  For  the  seagrass 
stations,  71%  were  correctly  identified  and  29%  were 
misidentified  as  bare  sediment.  Similarly,  83%  of  the 
sediment  pixels  were  correctly  identified  and  only 
17  %  were  misidentified  as  seagrass.  No  commission 
errors  occurred  for  algae.  For  algae,  50%  of  the  4  sta¬ 
tions  were  misidentified  as  either  sediment  (25%)  or 
seagrass  (25%).  However,  the  technique  correctly 
identified  algal  stations  along  the  central  Exumas,  the 
largest  algal-dominated  region  that  was  sampled. 
Applying  this  model  throughout  the  satellite  image 
(Fig.  7A),  seagrass  accounted  for  4.1  x  1010  gC  d_1  and 
covered  roughly  40  000  km2,  while  algae  accounted  for 
3.3  x  109  gC  d_I  and  covered  only  4500  km2.  A  major¬ 
ity  of  the  GBB  (50  000  km2)  was  comprised  of  bare 
sediment. 
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Fig.  7.  (A)  Mapped  distributions  of  algae  (either  macro-  or 
microalgae),  seagrass,  and  bare  sediment  on  the  Great 
Bahama  Bank.  (B)  Classification  accuracy  matrix:  diagonal 
values  (bold)  indicate  the  number  of  stations  that  were 
accurately  classified  in  each  category.  Values  outside  the 
diagonal  were  misclassified  stations  either  due  to  errors  in  the 
algorithm  or  variability  in  the  spatial  scale  between  single 
station  measurements  and  estimates  derived  over  the  1  km2 
pixel  area 

DISCUSSION 

Field  measurements  from  the  GBB  revealed  several 
different  types  of  benthic  habitats:  seagrass,  brown 
and  red  macroalgae,  and  sediment  with  and  without 
microalgal  biofilm.  Sediment  with  little  microalgal 
biofilm  or  seagrass  was  considered  the  prevalent  ben¬ 
thic  habitat,  covering  nearly  50  000  km2  of  the  bank. 
Seagrass  meadows  were  found  to  be  the  dominant  pro¬ 
ducer  on  the  bank.  Moderate  to  dense  stands,  primar¬ 
ily  of  turtle  grass,  were  mapped  over  37  000  km2  of  the 
bank  and  occurred  primarily  along  the  perimeters 
of  the  GBB  and  throughout  the  southern  bank.  Due 
to  high  winds  and  logistical  factors,  we  were  unable 
to  directly  sample  the  southern  portions  of  the  bank 
where  much  of  the  productive  seagrass  meadows  were 
identified  from  space;  future  field  campaigns  to  this 
remote  central/southern  section  of  the  bank  will  be 
required  to  validate  our  map. 

Using  the  lower  bound  estimate  of  growth,  NPP  esti¬ 
mated  for  the  densest  seagrass  meadows  was  on  the 
order  of  0.8  gC  m-2  d-1.  NPP  estimates  based  on  in  situ 


leaf-punch  studies  were  in  reasonable  agreement  with 
the  laboratory  measurements  of  leaf  metabolism  and 
varied  by  <20  %  for  different  ranges  of  LAI.  The  higher 
rates  of  NPP  measured  in  2002-2003  near  Lee  Stock¬ 
ing  Island  (see  Fig.  2B)  were  comparable  to  measure¬ 
ments  by  Koch  &  Madden  (2001)  of  1.65  to  2.29  gC  m-2 
d_1  in  a  lagoon  on  Grand  Bahama  Island  (Little  Bahama 
Bank)  using  in  situ  benthic  chambers. 

Benthic  algae  from  the  GBB  included  brown  and  red 
macroalgae,  as  well  as  microalgae  associated  with 
sediment.  Benthic  diatoms  have  been  identified  as 
the  dominant  microalgae  associated  with  all  sediment 
types  including  peloidal,  oolitic,  and  grapestone  sedi¬ 
ments  (Stephens  et  al.  2003).  Measured  rates  of  algal 
NPP  varied  widely  depending  on  the  type  and  concen¬ 
tration  of  algae  present  on  the  seafloor.  NPP  of  brown 
macroalgae  found  in  thick  3  to  10  cm  drifts  in  the  Exu- 
mas  regions  (Dierssen  et  al.  2009b)  was  very  high, 
upwards  of  18  to  30  gC  m-2  d_1  depending  on  macroal- 
gal  density.  In  contrast,  NPP  of  grapestone  with  rela¬ 
tively  thin  biofilm  was  barely  net  autotrophic  (0.005  gC 
m-2  d_1).  Average  NPP  for  tropical  macrophyte-domi¬ 
nated  ecosystems  has  been  reported  between  these  2 
extremes  at  -37  mol  C  m-2  yr_1  or  roughly  1.2  gC  m-2 
d_1  (Gattuso  et  al.  1998). 

Within  the  confines  of  this  multi-spectral  analysis,  no 
distinct  spectral  signature  was  found  to  effectively  dif¬ 
ferentiate  micro-  and  macroalgae.  Derivative  analysis 
showed  several  spectral  wavebands,  particularly 
between  505  and  550  nm,  where  various  producers 
could  be  differentiated  based  on  spectral  slope 
changes  that  may  be  useful  should  higher  spectral 
resolution  imagery  become  routinely  available.  With 
additional  spectral  information  (Phinn  et  al.  2008)  and 
other  coincident  measurements  and  analysis  (e.g. 
LIDAR,  texture  analysis,  etc.),  it  may  be  possible  to  fur¬ 
ther  refine  the  remote  sensing  tools  to  better  estimate 
the  vastly  different  productivity  rates  of  micro-  and 
macroalgae.  This  analysis  suggests  that  highly  produc¬ 
tive  algae  cover  roughly  4500  km2  of  the  seafloor  and 
represent  7  %  of  the  total  NPP,  but  this  number  may  be 
on  the  low  side  when  considering  ephemeral  blooms  of 
algae  and  with  inclusion  of  sediment  algae  within  the 
seagrass  meadows. 

Total  NPP  estimated  across  the  entire  GBB  was  5  x 
1010  gC  d_1  or  1.8  x  1013  gC  yr_1.  Propagated  error  is 
difficult  to  constrain  given  the  level  of  complexity  in 
the  input  data  and  model  formulation.  However, 
results  from  Fig.  6D  suggest  that  the  NPP  estimate  may 
be  accurate  to  roughly  within  a  factor  of  2.  Seagrass 
ecosystems  across  the  GBB  represent  around  0.04  % 
of  the  total  ocean  productivity  of  48  x  10 15  gC  yr_1 
(Longhurst  et  al.  1995).  If  these  rates  are  consistent 
across  the  -800  000  km2  of  shallow  water  carbonate 
sediment  banks  scattered  throughout  the  world  (Milli- 
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man  1993),  productivity  in  carbonate  systems  equals 
roughly  0.3  %  of  global  annual  ocean  productivity.  This 
is  consistent  with  previous  estimates  of  global  seagrass 
NPP  (Mateo  et  al.  2006).  Applying  these  techniques 
over  more  coastal  habitats  could  lead  to  better  con¬ 
strained  estimates  of  seagrass  to  global  NPP. 

Errors  in  the  NPP  estimates  derived  from  space- 
based  imagery  arise  from  many  sources  including 
atmospheric  correction,  bathymetry,  water  column 
optical  properties,  and  variability  in  estimating  NPP 
from  the  different  producers.  The  atmospherically  cor¬ 
rected  sea  spectral  reflectance  derived  over  the  bank 
was  representative  in  magnitude  and  spectral  shape 
across  a  variety  of  dark  and  bright  stations  sampled 
field  (see  Fig.  6B).  For  this  analysis,  bathymetry  was 
digitized  from  charts  because  good  agreement  could 
not  be  reached  between  in  situ  measurements  and 
optical  methods.  Bathymetry  maps  of  the  GBB  from 
SeaWiFS  (Stumpf  et  al.  2003)  and  Medium  Resolu¬ 
tion  Imaging  Spectrometer  (Lee  et  al.  2007)  imagery 
showed  many  of  the  same  features  reproduced  in  the 
present  study  and  these  techniques  hold  great  promise 
for  future  retrievals.  In  addition,  bathymetry  and 
benthic  properties  can  be  simultaneously  retrieved 
from  airborne  and  spaceborne  hyperspectral  sensors 
(Kutser  et  al.  2006,  Lesser  &  Mobley  2007). 

A  look-up  table  based  on  a  single  set  of  mean  water 
column  inherent  optical  properties  was  used  to  model 
bottom  reflectance  across  the  GBB.  Estimates  of 
bottom  reflectance  would  be  -15%  lower  and  NPP 
approximately  15  to  20%  higher  if  the  higher  attenua¬ 
tion  values  were  used  across  the  entire  bank.  The 
higher  attenuation  rates  west  of  Andros,  however, 
were  primarily  associated  with  fine  mud  sediments 
that  have  little  to  no  productivity  and  may  not  be  rep¬ 
resentative  of  the  water  column  overlying  benthos  with 
higher  NPP.  Dense  seagrass  meadows  can  reduce  cur¬ 
rents  inside  the  beds  (Scoffin  1970)  and  can  decrease 
the  particle  load  in  the  water  (Buonassissi  2009).  More¬ 
over,  the  inherent  optical  properties  used  in  the  radia¬ 
tive  transfer  modeling  were  in  the  range  of  those  mea¬ 
sured  around  Lee  Stocking  Island  (Boss  &  Zaneveld 
2003)  and  may  be  generally  representative  of  condi¬ 
tions  in  these  oligotrophic  waters,  although  some  vari¬ 
ability  in  absorption  by  colored  dissolved  matter  may 
be  evident  at  different  tidal  stages  (Otis  et  al.  2004).  A 
further  nuanced  approach  to  the  modeling  would 
incorporate  different  water  column  optical  properties 
expected  for  different  benthic  substrates  (Boss  &  Zane¬ 
veld  2003)  and  environmental  conditions  (e.g.  tides 
and  winds)  (Monahan  1993). 

In  addition  to  the  error  inherent  to  the  empirical 
regression,  estimates  of  NPP  may  also  be  in  error 
because  they  do  not  account  for  total  system-wide  NPP 
within  seagrass  meadows.  Seagrass  NPP  estimated 


from  leaf  growth  rates  represent  the  lower  bound  for 
seagrass  production,  because  they  do  not  include  the 
growth  of  below-ground  roots  and  rhizomes  or  pro¬ 
duction  by  other  autotrophs  associated  with  seagrass 
meadows.  For  these  types  of  seagrass  meadows, 
below-ground  productivity  may  be  as  much  as  35  %  of 
above-ground  production  (Duarte  &  Chiscano  1999). 

NPP  in  seagrass  systems  is  further  enhanced  by  the 
larger  autotrophic  community  consisting  of  phyto¬ 
plankton,  epiphytic  algae,  and  benthic  macro-  and 
microalgae.  In  these  clear  waters,  mean  phytoplankton 
biomass  estimated  from  the  2004  field  season  was 
0.12  mg  chi  a  m-3,  which  roughly  translates  to  primary 
production  of  9  gC  m-2  yr_1  (Behrenfeld  &  Falkowski 
1997),  an  insignificant  contribution.  Seagrass  system- 
wide  productivity  may  also  be  altered  by  epiphytic 
growth  on  the  seagrass  leaves  and  microalgae  within 
the  sediment.  Epiphytes  do  not  act  merely  as  neutral- 
density  filters  situated  above  the  seagrass  leaves,  but 
compete  for  photons  with  the  underlying  leaves  (Drake 
et  al.  2003,  Fyfe  2003).  Turtle  grass  epiphytes  from 
near  Lee  Stocking  Island,  Bahamas,  absorbed  a  maxi¬ 
mum  of  36%  of  incident  light  in  peak  chlorophyll 
absorption  bands  and  negatively  impacted  spectral 
photosynthesis  of  their  seagrass  hosts.  For  the  oxygen 
evolution  measurements  in  this  analysis,  however,  epi¬ 
phytes  were  left  intact  and  NPP  estimates  should  be 
inclusive  of  both  leaf  and  epiphyte  productivity  rates. 
Benthic  microalgae  from  the  upper  layer  of  sediment 
are  highly  variable  across  the  Bahamas,  as  noted  pre¬ 
viously,  and  some  productivity  is  expected  to  occur 
within  sediment  from  the  seagrass  meadows  them¬ 
selves.  Consequently,  system-wide  production  within 
seagrass  may  be  over  2-fold  greater  than  the  above¬ 
ground  seagrass  production  reported  here. 

Finally,  seagrass  photosynthesis  in  seawater  is  inher¬ 
ently  C02-limited  and  seagrass  productivity  may 
increase  in  response  to  rising  C02  levels  in  the  atmos¬ 
phere,  increasing  both  the  flux  of  organic  carbon  to 
the  sediments  and  the  rates  of  carbonate  dissolution 
(Zimmerman  et  al.  1997,  Invers  et  al.  2001,  Palacios  & 
Zimmerman  2007).  Rising  atmospheric  C02  may  initi¬ 
ate  a  positive  feedback,  further  enhancing  the  seques¬ 
tration  of  C02  into  bicarbonate  by  promoting  dissolu¬ 
tion  of  carbonate  sediments  on  shallow  banks  that  are 
heavily  vegetated  by  seagrasses  (Burdige  et  al.  2008). 
If  the  vast  majority  of  seagrass  biomass  decomposes  in 
the  sediments,  seagrass  productivity  across  the  GBB 
alone  may  sequester  as  much  as  2.4  x  1013  gC  yr_1  (sum 
of  above-  and  below-ground  productivity),  represent¬ 
ing  nearly  1  %  of  the  annual  oceanic  C02  uptake  of 
2.2  x  1015  gC  yr_1  (Feely  et  al.  2001).  If  these  rates  are 
consistent  across  all  carbonate  sediment  banks 
(-800  000  km2;  Milliman  1993),  seagrass  productivity 
may  be  responsible  for  as  much  as  8  %  of  the  ocean's 
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annual  net  C02  uptake.  Transient  macroalgal  blooms, 
on  the  other  hand,  may  be  associated  with  wind-driven 
processes  along  the  Exumas  (Dierssen  et  al.  2009b), 
and  their  abundance  has  been  associated  with  nutrient 
enrichment  in  coastal  estuaries  (Valiela  et  al.  1997). 
Benthic  microalgal  communities  serve  to  provide  food 
to  higher  trophic  levels  and  modify  nutrient  exchange 
between  the  water  column  and  sediments  (MacIntyre 
et  al.  1996,  Berelson  et  al.  1998).  Being  able  to  monitor 
the  present  and  future  distributions  and  productivity  of 
these  various  benthic  producers  will  be  important  as 
these  regions  are  influenced  by  both  natural  and 
anthropogenic  changes.  Enhanced  spatial  and  particu¬ 
larly  spectral  resolution  imagery  from  satellites  may 
afford  better  capabilities  for  such  future  monitoring 
and  analysis. 
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